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The spin excitations in the antiferromagnetic phase of half-filled Kondo lattice model are studied 
by means of the decoupling approximation for spin Green's function. The spin-wave spectrum 
is calculated as a function of Kondo coupling, and this is used to calculate the thermodynamic 
quantities at low temperatures. The Neel temperature of the form Tn = 0.087J 2 In (1/J) is 
obtained for the 3-dimensional case in the weak-coupling limit. It is pointed out that the ratio 
of the spin- wave velocity to the Neel temperature Vs/Tn is enhanced as the Kondo coupling 
becomes small, reflecting the long-range nature of effective interactions between localized spins. 



The electronic states of heavy-fermion compounds have been extensively studied to un- 
derstand a number of their anomalous properties due to strong correlations. The Kondo 
lattice model (KLM) and the periodic Anderson model (PAM) have been considered to give 
a microscopic basis to describe the characteristic features of these systems. It has been 
expected that KLM and PAM exhibit various types of magnetic long-range order due to 
effective RKKY interaction between localized spins. In particular, the properties of antifer- 
romagnetically ordered (AFO) phase realized in the half- filled system have been studied by 
the Gutzwiller variational method Q |Q and the slave-boson mean-field approximation |3| 
in connection with the magnetic instability of heavy-fermion semiconductors ||. However, 
the spin excitations and the resulting thermodynamic properties in AFO phase have not been 
fully understood, because the low-lying excitations are closely related to the non-local fluc- 
tuations neglected in these approximations. In the present paper, we are going to investigate 
these spin excitations in the AFO phase of KLM, making use of the spin Green's functions 
in the decoupling approximation similar to the Tyabrikov method for the Heisenberg model 
||. The spin- wave excitation spectrum and the Neel temperature are determined from these 
Green's functions. 

The KLM with the nearest-neighbor hopping is described by the Hamiltonian 

H = - t J2 4c jct + j]>>-4 (i) 

where c\ a and Cj CT represent the creation and annihilation operators of a conduction electron 
in Wannier state at the i-th site with spin a, respectively. The spin operator of conduction 
electron c?j is expressed in the form c?j = 1/2 C L 

Too'C-ia 1 , where t is the Pauli matrix. In 
the following, we set t = 1, measuring all energies in the unit of t. 

In order to describe AFO phase in this model, we introduce the staggered magnetizations 
of localized spin and conduction electron, respectively : 



m s for A sublattice 
-m s for B sublattice, 
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— m c for A sublattice 
m r for B sublattice. 



Then the Hamiltonian (1) is devived by the mean- field part Hq and the residual-interaction 
part H\ = H — Hq. The Hq describes the conduction band in the staggered field, and can be 
regarded as the quasi-particle bands with spin-density-wave (SDW) modulation by means of 
the canonical transformation 

4? + (2) 

with 



where c^ B ^ is the Bloch transform of cf} B ^ and = \J e^ 2 + (Jm s /2) 2 stands for the quasi- 
particle energy, and the wave vector A; runs over reduced Brillouin zone. Making use of such 
a transformation, the Hamiltonian (1) is rewritten in the form : 



where 



and 



H — Hq + H i , 

hq = £(^4^4? - ^4; )f 4; } ) - ^E(^ - ^) ( 4 ) 



fccT 



ffi = J^{(m, - S Az )(m c - a Az ) + (m s - Sf z )(m c - af z )} 

i 

+jJ2i< + S A q + <r}S% + a^ + S B - + o?-S*+}. ( 5 ) 

where o" g and 5 g are the Bloch representations of spin operators of conduction and localized 
spins, respectively. The simple mean-field treatment, which neglects Hi, produces the correct 
ground state of the large-5 limit, providing the appropriate basis for further analysis. In order 
to study the effects brought about by the fluctuation H\, we introduce the Green's function 
of localized spins, (S^ B ^ + ; S A q ), which satisfies the following equations of motion : 

z(st + ;S A -} = 2{s Az ) - JYl[(^ z st^;s A ;) + (S^ q ,a A , + ; S A q )}, (6 - a) 

q' 

z(S* + ; flf-> = - J£[<af Sf_+ ; Sf"> + (S B f q ,o B+ ; S A q )]. (6 - b) 

Then the decoupling approximations such as 

(a z S + ;S-)^(<r z )(S + ;S-), 

(S z a + ;S-)^(S z )(a + ;S-), 

are employed for the higher-order Green's functions, assuming the SDW state for carriers 
and Neel state for localized spins, respectively. This procedure is similar to the conventional 
Tyablikov method for the Heisenberg model. Thus we obtain 

(z - Jm c )(S A+ ; S A ;) = 2m s + Jm a Y^[44^(°™°w 
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+ U k U k+ q ( a t\ )ta k+q^ S -q >]> ( 7 - °) 

(z + Jm c )(5f + ; Sf"> = Jm s Y\ u k u k + M$ ] a i~l q v S A ~) 

k 

+« 9 <4! )t 4; ) 9i ;^->]- (7-6) 

In the above expression, the contributions from Green's functions of the type (a^aa^; S A ~) 
are neglected, because the effects of the particle-particle and hole-hole scattering are suffi- 
ciently small at the weak coupling and low temperatures. Employing similar approximations, 
the following equations are derived for the Green's functions appearing on the r.h.s. of eqs.(7- 
a,b) 

(z±E k ±E k+q )(aifa^ ql ;S A q } = 

^{1 - f(E k ) - f(E k+q )}[u± u ± +q (S A+ ; Sf"> + uM +q (S* + ; ( 8 ) 

where f(E k ) is the Fermi distribution function. From this closed set of equations, we obtain 
the explicit expression for the spin Green's function 

(S q A+ ;S A ;) = 



2m, 



{1 - ^ ff (*)}z - ^{X AA (0) + X AB (0) ~ X AA (z)} 



{1 - ^ q (z)} ^ - (^) [{x AA (o) + x AB (o) - X AA (z)V - {x AB (z)Y] 



(9) 



with 



where use has been made of the relation 

m c = -^{ X AA (0) + X AB (0)}- (11) 

Replacing z in eq.(9) by u + i8, the selfconsistent equation for the staggered magnetization 
m s is given by 

]T/ dcnHlm(^ + ;^-), (12) 



where n(w) is the Bose distribution function (e^ — l)^ 1 and j3 = T _1 . 

This Green's function (9) has two types of poles corresponding to the Stoner and spin- 
wave excitations, respectively. Since the energy band of spin-wave excitations characterized 
by the energy scale 2Jm c always lies lower than the continuum of the Stoner excitations 
E k+q + E k which have a gap 2Jm s in the whole momentum space, in the weak coupling 
regime the spin-wave dispersion is well approximated in the form 



flm, VW" 4 + xi B - xi A Y - {xf 3 } 2 



A <*> = V r-jw • (13) 
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by replacing \ (z) and 4> ( z ) by their values at z = 0. Our result (13) clearly shows the presence 
of the infinitesimal spin excitation at q=0 which is expected from the Goldstone theorem. 
We note that in the single-site approximation such as Gutzwiller and slave-boson methods, 
there remains a finite spin excitation gap even in AFO phase. Since the spectral intensities 
of Stoner excitations are sufficiently small compared to those of spin-wave excitations in the 
weak couling regime, in the following analysis we neglect the contributions from the Stoner 
excitations. 

Let us first consider the asymptotic form of this dispersion at zero-temperature in the 
weak-coupling limit J — > 0. In this limit, we can show that Xq A an d Xq B diverge as In J, 
while Xq A an d X AB (q 7^ 0) tend to a constant value. Therefore the equation (13) becomes 
independent of q in the limit J — > with a jump at q = 0. This means that the results 
of the present theory at zero-temperature reduce to those of the mean-field theory in the 
weak-coupling limit, even although the effects of quantum fluctuation are included in the 
present decoupling scheme. This asymptotic behavior is independent of the dimensionality 
of the system and consistent with our previous work of ID PAM j?J . 

For general couplings, equation (13) has been solved numerically, and the results are 
shown in Fig.l: Here, considering the weak coupling case, we use m s = 1/2 and neglect devi- 
ation from it. The spin-wave dispersion strongly depends on the strength of Kondo coupling. 
It is found that the rapid increase in small-q region is gradually relaxed as J becomes large. 
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Fig.l. Spin-wave dispersion, A (q) / A(n/2), as a function of q = (q,q,q) for various Kondo 

coupling. The broken line represents the result in Heisenberg antiferromagnet ; yl — (cos (q x ) + cos (q y ) + cos (g 2 )) 2 /9. 
The linear expansion in q for eq.(13) yields the spin- wave velocity 



v s = 0.088jyin^-J, (14) 

for 3D KLM in the weak-coupling region. Similar results are also obtained for the lower 
dimensional systems : v s /Jy/ln (l/J) = 0.135 for 2D and 0.138 for ID. Once the spin-wave 
velocity is determined, we can analytically calculate the specific heat and parallel suscepti- 
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bility in the low-temperature limit || || as follows 



XII 



4V3tt (T 



(15) 



(16) 



We can obtain the Neel temperature Tn by linearizing eq.(12) with respect to m s . Thus 
the equation for Tn is given by 



Tn 
J 2 



Xn(o) 



2 2 Xn (0) + Xn(9) 
^ a Xn(0) - XN(ff) 



n -l 



where 



Xn(?) 



2 

iV 



E 



1 - /N(efc) - fl<s(£k+q) 



(17) 



(18) 



and /N(efc) is the Fermi distribution function (e" £ft + 1) at T = Tj$. Equation (17) has 
been solved numerically with respect to Tn for various Kondo couplings J, and the results 
are shown in Fig. 2. It is found that the Neel temperature shows a monotonic increase with 
increasing J. 




Fig. 2. Neel temperature Tn as a function of Kondo coupling. The solid and broken lines 
represent the present and the mean- field results, respectively. The inset shows that the Neel 
temperature is linear to J 2 In (1/J) in the weak-coupling region. 

From the asymptotic behavior shown in the inset of Fig. 2, we obtain the expression of Tn 
in the weak-coupling limit as follows : 

T N = 0.087 J 2 In (jj . (19) 

The mean-field theory also predicts the J-dependence similar to (19) in the limit J — ► 0. 
However Tn in our theory is smaller than that in the mean-field theory by a factor 0.71. The 
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deviation from 1 is due to the quantum fluctuations between localized spins. This is to be 
compared with the corresponding ratio 0.66 in Heisenberg antiferromagnet. 

It should be pointed out that the ratio of the spin-wave velocity and the Neel temperature 
v s /Tn is enhanced as the Kondo coupling becomes small. This is because the long-range 
nature of effective interactions between localized spins suppresses collective excitations at 
low temperatures. 

In summary, the spin excitations of the AFO phase in half-filled KLM were studied by 
means of the decoupling approximation for spin Green's function. The spin-wave spectrum 
and the Neel temperature were calculated as a function of Kondo coupling, and their asymp- 
totic behaviors in the weak-coupling limit were discussed in detail. It was shown that the 
collective spin excitations in weak-coupling KLM are considerably suppressed in comparison 
with those in the Heisenberg antiferromagnet, reflecting the long-range nature of effective 
interactions between localized spins. 
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